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Abstract The dynamical coupling between the solar chromospheric plasma and 
magnetic field is investigated by numerically solving a fully self-consistent, two- 
dimensional initial-value problem for the nonlinear collisional MHD equations in¬ 
cluding electric resistivity, thermal conduction, and, in some cases, gas-dynamic 
viscosity. The processes in the contact zone between two horizontal magnetic 
helds of opposite polarities are considered. The plasma is assumed to be initially 
motionless and having a temperature of 50,000 K uniform throughout the plasma 
volume; the characteristic magnetic field corresponds to a plasma /3 > 1. In a 
physical-time interval of 17 seconds typically covered by a computational run, the 
plasma temperature gradually increases by a factor of two to three. Against this 
background, an impulsive (in 0.1 seconds or less) increase in the current-aligned 
plasma velocity occurs at the site of the current-layer thinning (sausage-type de¬ 
formation, or m = 0 pinch instability). Such a “velocity burst” can be interpreted 
physically as an event of suprathemal-proton generation. Further development 
of the sausage instability results in an increase in the kinetic temperature of the 
protons to high values, even to those observed in flares. The form of our system of 
MHD equations indicates that such increases are a property of the exact solution 
of the system at an appropriate choice of the parameters. Magnetic reconnection 
does not manifest itself in this solution: it would generate flows forbidden by the 
chosen geometry. Therefore, the pinch-sausage effect can act as an energiser 
of the upper chromosphere and be an alternative to the magnetic-reconnection 
process as the producer of flares. 
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1. Introduction 


The chromosphere is much less studied than the underlying photosphere and 
overlying corona. However, its properties are widely discu ssed i n the context of 
the f undamental problems of solar physics lAschwanden. [20o3 : Fletcher et ai, 
l2010l ). Much effort has been made to understand the mechanism of heating of the 
upper chromosphere and corona. Such a mechanism is usually associated with 
coronal processes at plasma-parameter values /3 < 1. On the other hand, some 
observed phenomena lead, as Aschwanden notes, “to a paradigm shift from the 
coronal he ating p roblem to a dynamic chromospheric energisation problem” (As¬ 
chwanden, I 2 OOI b the latter cor respon ds to the case of /3 > 1 (Aschwanden and 
Nightingale, 20051 : Aschwanden. |200^ . Therefore, particular attention should be 
given to plasma properties at /3 ~ 1. 

Progress in understanding these issues will remain slow as long as the knowl¬ 
edge of the chromosphere remains fragmentary. In many cases, detailed in¬ 
formation provided by the modern high-resolution instruments can hardly be 
compared with theoretical conclusions, since the theories are based on non¬ 
observable parameters (nonpoten tial m agnetic field, magnetic stress or twist, 
electric fields, etc.) 1 Aschwanden. l200^ . 

Such difficulties in studying the dynamics of the chromospheric plasma itself 
could be avoided using the completely self-consistent system of MHD equations. 
Various quantities are mutually related in a solution of this system, and a passage 
from some quantity crucial for the process to another, more convenient for an 
observer, can be done quite routinely. 

However, in the practice of simulations of the chromospheric plasma, even the 
incl usion of magnetic fields increases enormously the level of complexity (Carls- 
son, 2007t) . Highly dynamic, nonlinear processes are typical of the chromosphere. 
Small scales are involved in the m due to both turbulent chaos and the formation 
of shock fronts. Aschwanden ( 2004[ ) writes that “typical difficidties with the 
MHD method are the same as they are typical for numerical simulations, such 
as heavy computing demand, convergence problems, insufficient spatial resolu¬ 
tion to handle discontinuities, line-tying (continuous slippage of magnetic-field 
lines)”. 

However, direct simulations of dynamical coupling between laboratory ion- 
electron plasma and magnetic field were made by numerically solving a two- 
dimensional initial-value problem for fully self-consistent, nonlinear collisional 
MHD equations, which give a un ified descripti on for /3 > 1, /3 Ri 1, and /3 < 1 
(Brushlinskii and Morozov, 198fll : Brushlinskv. IlfiSOli . Later, various situations 
were considered in the cont ext of studying laboratory plasma-channel flows (see 
an overview by 1 Alekseeva. l2006h b It is worthwhile to try to comprehend how 
the features of the dynamical coevolution of the plasma and magnetic field, 
revealed in this way, can manifest themselves under solar conditions. To this 
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end, we solve here numerically an initia. l-value problem of this sort; in so doing, 
as in Brushlinskii and Morozov ( 1980ll . we do not assume any properties of 
the solutions sought for. In view of investigating the solar-plasma dynamics, we 
employ a substantially improved numerical technique. This makes it possible to 
describe small-scale turbulence and discontinuities, automatically passing, when¬ 
ever necessary, to the use of the class of generalised functio ns. Th is technique 
was v erified in our previous studi es (Kshevet skii and Gavrilov. l2005l : Kshevetskii, 
20061 : Alekseeva and Kshevetskii, 20061 . 201lh 'l. This allows us to obtain solutions 


for comparatively long time intervals even if small-scale perturbations develop. 

We consider a two-dimensional problem assuming that all physical quantities 
are constant along the magnetic-field lines and the magnetic field is horizontal. 
This is motivated by the fact that about 95 % of the magne tic flu x issuing from 
the p hotosphere closes below the coronal heights (Priest, Il982l : Aschwanden, 
20041 1. so that the magnetic field in the chromosphere consists largely of hori¬ 


zontal apical segments of the field lines; horizontal segments are also present in 
the field lines of the magnetic canopies. 

We investigate here the coevolution of the collisional electron-proton plasma 
and magnetic field at the characteristic plasma parameter val ues = 1-5, 
1.6, or 2.3, which are presei it in t he upper chromosphere fGary. [2001 1 [see also 
Figure 1.22 in Aschwanden (120041 1]. Our MHD simulations are carried out under 
the assumption that the plasma is initially motionless and has a temperature of 
50,000 K everywhere in the computation domain. 

Such fully ionized plasma can be attributed to the upper chromosphere. 
The upper boundary of the chromospheric layer is typically associated with a 
temperature of about 50,000 K; this layer is assumed to be overlaid with a tran¬ 
sition region, which, in tu rn, changes into the corona at temperatures of about 
500,000 K IGabriel. Il994 1. However, more recent observations modified our un¬ 
derstanding of the geometry of the solar atmosphere, which is now regarded 
as an inhomogeneous mixing of photospheric, chro mosp heric, and coronal zones 
produced by variou s dynamic processes fSchriiver. l200lll [see also Figure 1.17 of 
Aschwanden ( 2004[ l]. Since we consider small-size regions, there is generally no 
need to attribute them to particular layers, and the temperature can be assumed 
to be the basic parameter cha racter ising the plasma under studyQ 

Alekseeva and Kshevetskii ( 201lll carried out simulations for a characteristic 
plasma- /3o* = 1.5 and an initial magnetic configuration consisting of two re¬ 
gions with oppositely directed magnetic fields. We noted the development of a 
pinch-type plasma instability. Some its features are also observed in laboratory 
plasmas; these are the contraction of the current layer as a whole, the local 
thinning (sausage-type deformation) of the stagnated current layer, and the 
development of oppositely directed, current-aligned plasma streams. However, 
the process in an unbounded solar plasma proved to be more complex. An¬ 
other pair of jets accompanied by a strip of small-scale inhomogeneities was 


^To avoid misund erstan ding, it is worth noting that the definition of the upper chromosphere 
used by Gabriel ll 19941) is not the one in common usage, and there are also other different 
definitions, usually based on the characteristics of neutral-particle distributions over the solar 
atmosphere. We can generally consider plasma with temperatures of above 20,000 K using the 
equations for the fully ionised gas. 
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directed transversely to the electrical current. As viewed along the magnetic- 
held lines, the pattern of field, velocity, and density distributions appeared to be 
cross-shaped; we shall call such a pattern a transient “cross-shaped” strueture. 
The jet velocity exceeded 20 kms“^. A class of situations was identihed in 
which the plasma temperature increased in the central part of this pattern (the 
computations were terminated when it doubled). 

Here, we shall demonstrate that impulsive processes can develop against the 
background of a gradual development of this pattern (which will be illustrated 
in the case of /3o* = 1-6). Specifically, intense plasma streams emerge in a 
temporal interval of about 0.1 second. Such a velocity burst can physically be 
interpreted as an event of suprathermal-ion generation. The kinetic temperature 
based on the maximum velocity can reach about 0.5 MK. This is indicative 
for the dynamic energisation of a parcel of chromospheric plasma to coronal 
temperature. We conclude that the velocity burst is due to a sausage-type pinch 
instability. This velocity peak is neither unique nor the highest one. The very 
form of the initial-value problem considered here indicates an important property 
of its exact solution: the development of the sausage instability ultimately results 
in an infinite growth of the plasma velocity [?;] with time [t], which implies an 
increase in the kinetic temperature to high values, even such as observed in 
flares. Our numerical solutions allow us to trace this process. 

The infinite velocity growth can naturally be suppressed by gas-dynamic 
molecular viscosity, which will be observed in simulations with a strongly ex¬ 
aggerated viscosity. Some other factors may also act similarly. Analyses of their 
role fall beyond the scope of our study. We are interested in investigating the 
possibility of impu lsive plasma energisation by the pinch instability and, follow¬ 
ing Carlsson ( 2007l b “perform numerical experiments and modelling in simplified 
cases in order to fashion a basic physical foundation upon which to build our 
understanding.” 


2. Magnetohydrodynamic Background 

2.1. Collisional Plasma and its Dynamical Properties Depending on 
the Rarefaction Degree 


Some particular phenomena known in the theory of accelerating plasma channels 
may be of interest to solar physics. Two-dimensional flows transverse to the 
magnetic field with the physical quantities constant along the magnetic-field lines 
have long been studied. Early numerical simulations ~ generally, with a finite 
electric conductivity and the Hall effect taken into account -discovered the 
development of closed solitary electr ic-current structures, i.e. isolated ma gnetic 
tubes (Alekseeva and Solov’ ev, Il964 : Brushlinskii, Gerl akh, a nd Morozov, 1968 : 
Morozov and Solov’ev. Il98fll : Brushlinskii and Morozov, 198(ll l. Without the Hall 
effect, they develop at /3 > 1, and they do n ot dep end on the presence of walls 
confining the plasma (Alek seeva and Solov’ev, 1964 : Morozov and Solov’ev. ll98Cll : 
Brushlinskii and Morozov. Il98fll l. The role of the Hall effect in the dynamics of a 
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dense plasma with a magnetic field is determined by the dimensionless parameter 

^ = c(eL)“^\/TOi/(47riV*), (1) 

which will be referred to as the Hall-plasma-dynamics parameter; here, c, mp 
L, and are, respectively, the speed of light, the proton mass, the length 
scale, and the reference concentration of particles of the a given electrical-charge 
sign. The plasma is assum ed to be quasi-neutral, with a large but finite 
(Brushlinskii and Morozov, 1980ll . In a Hall plasma, although structures of the 
above-mentioned sort arise at /3 > 1, they manifest themselves more clearly 
(and with sharper gradients), as /3 decreases, being rela ted to explosive events 
characteristic of Hall plasmas (Brushlinskii and Morozov, 1980ll . These phenom¬ 
ena are, to all appearance, due to a gl obal g radient of the gas pressure in the 
flow direct ion (A le kseeva and Solov’ev. Il964 : Brushlinskii and Morozov, 1980l : 
Alekseeva, 1997a 1999ll . 

In the solar atmosphere, the Boltzmann-distribution-based pressure gradient 
can produce similar phenomena. A stationary two-dimensional solution of the 
same equation system was found analytically for the case of a very w eak magnetic 
held in a Hall plasma with the presence of gravity (Alekseeva, 200fil l. It represents 
a thin layer of a current directed upward, with a downward plasma how and a 
slower return how at the periphery of this layer. It changes into a thin current 
sheet in a more rarehed plasma, which may be of interest in the context of 
magnetic reconnection. Goodman ( 2005[l used a steady-state MHD model of 
the collisional chromospheric Hall plasma, with the electrical-conductivity and 
thermolectric tensors taken into account in Ohm’s law. For special conditions of 
a given magnetic held of a current-layer type and a given bulk plasma how 
orthogonal to the magnetic held, he found that the plasma thermodynamic 
parameters in a certain part of the current layer correspond to the lower-corona 
range. 

These results prompted us to start a systematic numerical investigation based 
on the collisional MHD equations for fully ionised chromospheric layers, generally 
with the inclusion of gravity, the Hall effect, and time-dependent electric and 
thermal conduction. To solve an axisymmetric or planar-symmetric initial-value 
probl em, we applied a special numerical method (Kshevetskii and Gavrilov, 
20051 1 that uses generalised functions wherever necessary, to the solar plasma. 


2.2. MHD Chromospheric Processes in our T-wo-Dimensional 
Numerical Simulations 


The penetration of the magnetic held into the chromosphere immediately and 
typically entails i) a tendency of the gas pressure and magnetic pressure to 
come into balance and ii) the interaction between the penetrating magnetic 
held and other closely located chromospheric magnetic helds. Generally, the two 
phenomena parallel each other and both produce wavelike steadying processes 
(from here on, referred to as primary waves), which interfere. 

We investigate this typical situation for the case of two contacting magnetic 
regions of opposite polarities. (For simplicity, we choose the initial distributions 
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of the variables to be symmetric with respect to the zero-field surface. We also as¬ 
sume the physical quantities to be constant along the magnetic-field lines, which 
are horizontal, straight and parallel; therefore, a 2D problem can be formulated.) 
In this circumstances, the primary waves associated above with phenomenon ii) 
are due to electrical-current contraction in the plasma medium and are formed 
by the compression and subsequent expansion of the whole current zone between 
the regions. Also, a sausage instability similar to that in the laboratory Z-pinch 
can be expected at some locations in the same zone, in the presence of the 
primary steadying waves and/or after they recede. 

Our numerical method allows us to calculate inhomogeneities (waves, turbu¬ 
lence, and shocks) spontaneously arising during the nonlinear process of plasma 
and magnetic-field coevolution. We do this without anticipating any properties of 
the solution. Thus, solving an initial-value problem for the fully consistent system 
of MHD equations, we can perform direct numerical simulations specifying only 
the initial conditions. 

Our early simulations at ^ 0 or, in other words, taking into account the 

Hall effect and gradient of electron pressure showed that t he pro cess modeled is 
highly dynamical and complex (Alekseeva and Kshevetskii, 2006ll and it is virtu¬ 
ally impossible to distinguish between the contributions of different phenomena 
to the observed scenarios. The difficulties can be overcome only by a stepwise 
consideration of progressively more complex situations. Here, as in Alekseeva 
and Kshevetskii ( 201ll l. we omit the Hall effect and electron-pressure terms in 
the equations [i.e. we use the system of equa tions referred to as the standard 
MHD model (see, e.g., Leake and Arber, 2006l l] . 

Thus, at least three groups of physical phenomena can manifest themselves in 
our direct simulations: primary waves of steadying processes, pinch sausage in¬ 
stability and turbulence of various kinds, including small-scale waves and shocks 
due to the nonlinearity of the processes. 

Simulations by Alekseeva and Kshevetskii (120 llh revealed a time interval 
during which cross-shaped transient structures are present. Here, we use a higher- 
accuracy algorithm to demonstrate that the transient structures can be accom¬ 
panied by impulsive generation of suprathermal protons in a manner similar to 
the formation of proton beams in a laboratory Z-pinch. 

To conclude this section, it is worth noting that our approach differs from 
those typically adopted. Many researchers avoid considering the effects of spon¬ 
taneously developing small-scale turbulence, waves and shocks because of se¬ 
vere difficulties entailed by such effortsU In contrast, our method is specifically 
tailored for describing small-scale inhomogeneities irrespective of their scale. 


^For example, Leake and Arber 11200(^1 . noting the difficulties with the inclusion of small-scale 
effects in direct large-scale MHD simulations, artificially introduce a sort of relaxation - in 
particular, to suppress unresolved acoustic modes, thermal conduction and shock dissipation 
on small scales. In this way, they succeed in their simulations with chromospheric neutral 
particles taken into account. 
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3. Mathematical Formulation of the Problem 


3.1. The System of MHD Equations 


To nondimensionalise the variables, we use two basic units, viz. the number 
density [N^] of either component of the quasi-neutral proton-electron plasma and 
the plasma temperature [T*], which is characteristic of a certain layer of the solar 
atmosphere before the arrival of the magnetic field and refers to the nonmagnetic 
regions after the field arrival. Then p* = mjiV* will be the unit density. We 
also choose the scale height of the atmosphere [H = kT^ / (m-ip) (where k is the 
Boltzmann constant and g is the gravitational acceleration)] as the unit length 
L. 

Our problem consists in studying the processes that can occur in such a 
layer under the influence of a magnetic field. As a unit magnetic field, we 
use its characteristic magnitude [Bo]- Then the combinations Bo^/(47r) and 
= Bo/-\/47rp» can be used as the unit pressure and velocity, respectively. We 
choose to* = as the unit time and cBq/(AttH) as the unit of the electric 

current. Thus, all of the units are specified, given the dimensional quantities T,, 
N^, and Bq. 

However, we should keep in mind that the response of the plasma medium 
to the initially present magnetic field (with a characteristic value of Bq) can 
conveniently be judged by the value of the plasma-,S parameter, 

/3o* = SttP./Bo^ (2) 


it gives a general characterisation of the situation, being the ratio of the char¬ 
acteristic dimensional pressures - the gas pressure P* = kN^T^^ to the magnetic 
pressure Bq/Stt. We write Bq and, accordingly, uq* in terms of /3o* (recall that 
T* and iV* do not depend on the magnetic field) as 

Bo = (3) 

z;o* = /3o*”^'^^\/2fcT*/mi, (4) 


to express all the introduced units in terms of T*, W, and /3o* • 

We restrict ourselves to a two-dimensional geometry and assume that the 
physical quantities do not vary along the straight, parallel, horizontal magnetic- 
field lines, and the plasma motion is directed across them. In this case, the 
complete system of self-consistent nonlinear equations of collisional magnetogas- 
dynamics with allowances for finite electrical resistance, ther mal c onductivity, 
and p lasma gas-dynamical viscosity (Brushlinskii and Morozov, 19801 : Braginskii, 


ipi 

19651 1. written in terms of the dimensionless functions B, T, p, V, v and the 


dimensionless electric current [j], has the form 


fdv 

\dt 


+ (^,.V)^; =-V B 


2 


MAv, 


f+ V.(pn) = 0, 


(5) 

( 6 ) 
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^ /^O* rri 

'P=—PT, 

(7) 

3B 

— = Vx(uxJB)-Vx (0j), j = V X B] 

(8) 



V- {KVT) + Qf. 

(9) 

5/3; the functions 


e = 

(10) 

K = 

(11) 

M = 

(12) 


are, respectively, the local magnetic diffusivity, dimensionless local thermal con¬ 
ductivity and dimensionless kinematic viscosity coefficient (assumed to be inde¬ 
pendent of the local temperature variations); the dimensionless factors 0*, k*, 
/r* do not depend on /3o*, being determined in Appendix A by the originally 
specified dimensional parameters of the medium - specifically, 0* depends on T,, 
and both k* and /i, depend on T, and N^. We do not take into account 
thermoelectricity and neglect the influence of gas-dynamic viscosity on plasma 
heating. As already noted, the term taking into account the effect of gravity is 
omitted in Equation (I^. The equation V • B = 0 is automatically satisfied for 
the assumed two-dimensional geometry. 

We will consider here the dependence of the variables on two spatial co¬ 
ordinates in a plane perpendicular to B, in which we introduce a Cartesian 
coordinate system (a;, z) with the x-axis directed vertically downward. The dis¬ 
tribution of B{x, z) over the plane is displayed as a greyscale map; we note that, 
in the geometry considered, the contours of the magnetic-field strength coincide 
with the electric-current lines. 

3.2. Initial Conditions 

Let the absolute value of the initial magnetic field be symmetric with respect 
to two symmetry axes, x = Xc = Xniax/‘2 [where B(xc, z)\t=o = 0] and 2 = 


®In their Sections 2.3 and 2.4, Brushlinskii and Morozov lll98Clf) use the system of Equa¬ 
tions 0-lll to investigate the plasma and magnetic-field dynamics without taking into 
account the Hall effect. Thus, this system describes an isotropically conducting plasma. If, 
however, this effect is included, the system should be supplemented with terms with a pa¬ 
rameter ^ defined by Equation QJ and characterising the difference between the macroscopic 
(dimensionless) velocities of t he ion and electron gases, Vi — Ve = /p; see Sections 1.2 and 2.5 

of Brushlinskii and Morozov (llQSOT) . The presence of these terms implies a tensor form of the 
transport coefficients, which is equivalent to the appearance of the ^-dependent Hall and 
Pedersen con ductiv ities, cth and cjp. Some important phenomena discovered by Brushlinskii 
and Morozov lll980fi could be accounted for in terms of the local variations of cth an d crp in the 
cours e of the plasma and magnetic-field coevolution (Morozov and Solov’ev, 119801 : Alekseeva, 
Il980l^ . 
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Zc = Zmax/2. The distribution of B(x, z)\t=o is shown in Figure [TJi, where the 
following notation is used: 

X = {x :rc)/Xniax; Z = (^Z Zq) / (T^) 

An analytic formula that specifies the initial magnetic field [B(x, 2)|t=o] shown in 
Figure[T^ is given in Appendix B. The initial field thus chosen is the sum of two 
terms (see Appendix B): the first term is smaller in absolute magnitude and is 
distributed over the light grey and dark grey areas in Figure[lji; the second one is 
greater and more localised (narrow black and white strips in the figure). Since the 
equation V • B = 0 is automatically satisfied in the two-dimensional geometry 
of the problem, each of these terms can be considered a separate magnetic field, 
either a distributed or a more localised one. 

We assume that the plasma is initially motionless and isothermal: 


v{x,z)\t=o =0, r(x, z)|t=o = 1- 


(14) 


In this case, according to Equation 0, the density at t = 0 differs from the 
pressure only by a constant factor. The initial condition for the pressure should 
specify, however, the degree of balance between the initial distributi ons o f the 
magnetic and gas pressure. Here, as in Alekseeva and Kshevetskii (120111 ). we 
restrict ourselves to two cases. 

The case of no initial balance (NIB). If the pressure of the initially specified 
magnetic field B[x,z)\t=o is not balanced by the gas pressure, then, according 
to our choice of units, 


p{x,z)\t=o = 1 

and, accordingly, in view of the Equation © of state, 


(15) 


V{x,z)\t=o = ^o*/2. 


(16) 


The case of partial initial balance (PIB). We assume that, in the case of 
PIB, the distributed magnetic field, being in balance with the local gas pressure 
P(a;,2)|t=o, is superposed with an unbalanced localised magnetic field and the 
balance thus appears to be only partial. Physically, this case can be interpreted 
as the penetration of a new, unbalanced magnetic field to the area of the old, pre¬ 
existing field (which has already reached a pressure balance). The corresponding 
formulas for the initial pressure and density are given in Appendix C. 


4. The Properties of the Exact Solution Dictated by the 
Formulation of our Initial-Value Problem. 

Our analyses of the form of Equations © - ® with initial conditions given by 
Equations (O, (EH) together with either Equations (fT5|) . (EH) or (EH, dsni), lead 
to the following conclusions. 

i) If the initial magnetic and gas pressures are completely balanced, a station¬ 
ary solution exists. However, in both the cases of a partial balance (PIB) and of 
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Figure 1. Distribution of B{x, z) in the plane (x, z) for the outset stage of the magnetic field 
- plasma coevolution at /3o* = 1.5 in the absence of gas-dynamic viscosity at different times: 
(a) t = 0, the cases of both PIB and NIB (see below in this section); (b) t = 0.008 and (c) 
t = 0.024, the case of PIB. The x-axis is directed downward from the upper left corner, and 
the z-axis is directed to the right. The quantities X and Z are related to x and z according to 
Equations la. The image is stretched vertically by a factor of 4.4. Light: positive magnetic 
field (directed inward); dark: negative magnetic field. The white dot on the dark background 
and the black dot on the light background mark the positions of the minimum and maximum 
B{x, z)\t—o, respectively. [The other greyscale maps represent the central part of this domain, 
marked with the white frame in panel (a).] 
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no initial balance (NIB), dynamical processes manifest themselves. We discuss 
them now assuming that M = 0 in Equation 

ii) If the initial configuration of the magnetic field has the form of Equa¬ 
tion (l25l) . illustrated in Figure [T^, the subsequent evolution of the magnetic 
field and plasma state should be similar to those in a laboratory Z-pinch. Note, 
however, that the contraction of the electric current to the axis x = Xc does not 
result here in substantial increases in the maximum \B\ in view of the planar 
geometry. 

iii) In the NIB case, the coevolution of the plasma and B begins with the 
origin of a system of powerful waves, fast electric-current contraction, collapse 
of the plasma and its subsequent retreat. All of these processes strongly affect 
the bulk of plasma. 

iv) In the PIB case, the early evolution is more quiet. The unbalanced magnetic-] 
field component d present between the layers of the balanced field (Figure [T]i) 
form “tongues” stretching to x = Xc, which provoke the sausage instability. 
Since B = 0 at a; = Sc, an a; component of the magnetic-pressure gradient 

is present near this axis; according to Equation (O, oppositely directed flows 
with a velocity Vx^ directed toward this axis from the both sides, develop. The 
form of Equation (jS]) indicates that the magnetic field is partially frozen in 
the plasma. Therefore, these flows carry the magnetic field closer to the axis, 
which results in a further steepening of the magnetic-field gradients, and so they 
intensify themselves. The substance can freely flow out along the neutral line 
and then spread from it (no electrodes or walls are present, in contrast to the 
conditions of laboratory setups). As to magnetic viscosity [0], it is not effective 
damper of the process, since the great gradients of B correspond to great j, i.e. 
enhanced heating of plasma and decreasing of 0 in accordance with Equations 
©mi). Therefore, in the absence of gas viscosity M, the initial-value problem 
of Equations (O-®, (HH), (1^ . (fTSl) . (ITCl) . or (OHl 15(11) has solutions with the 
velocity [u] approaching oo with t. The velocity increases the faster, the closer 
the magnetic “tongues” to the axis a; = Sc. It is difficult to predict whether this 
process will be monotonic or impulsive. 

v) The most important feature of the process considered is the presence of 
opposite plasma flows approaching each other through the zone where B k. Q. 

It is therefore reasonable to conjecture that the inclusion of the gas viscosity 
\M] in Equation ([5]) can slow down or even stop the process. (In contrast to the 
magnetic viscosity [0], the gas viscosity increases with temperature according to 
Equation (1231) 1. Here, we are interested in a qualitative evaluation of the role of 
the gas viscosity and ignore, for simplicity, its local variations with temperature. 


5. The Procedure of Numerical Solution 

The verification criterion for a numerical code solving the two-dimensional initial-] 
value problem in Equations (O-®, (HU), (051) . (flSl) . (fTO . or (1^ 15(11) is the 
similarity between the numerical and exact solutions. However, the numerical 
solution possesses some features that are not present in the exact solution. For 
example, the exact solution of the problem in Equations ®-®, (HU), (EH), (HH, 
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(HH), or dMiisni) must be of the type of symmetry relative to the axis x = Xc 
considered and cannot describe the wriggle (or kink) instability of the Z-pinch, 
whereas its numerical solution can; the reason is that numerical noise prevents 
accurately describing the symmetry imposed by the initial conditions. 

To solve our initial-value problem numerically, we use a finite-difference scheme| 
that approximates the equations with the second order of accuracy in time and 
space; central differences are used to approximate the spatial derivatives. It is 
known that stable numerical methods of second-order accuracy introduce no 
spurious dissipation (although they may introduce spurious dispersion; methods 
of a Hrst-order accura cy int roduce spurious dissipation without spurious disper¬ 
sion - see LeVeque, Il992l l: on this basis, many numerical algo rithms free of 
spurious dissipation have been developed (Thomas and Roe. 119931 ). Our second- 
order-accuracy technique belongs to this class of methods free of numerical 
dissipation. 

Our scheme is constructed so as to ensure the existence of a non-increasing 
grid analogue of the wave-energy functional; this immediately proves a stability 
theorem for the method. The scheme is conservative, i.e. the energy, mass and 
momentum are conserved. As already noted, the solution of the problem can 
be sought in the class of generalised functions; the algorithm is able to auto¬ 
matically make transitions from smooth to generalised solutions and vice versa. 
Therefore, the method makes it possible to obtain non-differentiable solutions, in 
particular, with discontinuities (even multiple). This property of our numerical 
method proves to be important, since we will study the breakdown of the initially 
smooth field into small-scale waves and other structures, i.e. into the features 
that can be described either by functions with large derivatives or even by non¬ 
smooth functions. Our numerical method generalises the previously developed 
method for the numerical integr ation of equations of motion of the atmospheric 
gas (Kshevetskii and Gavrilov, 200^. Structurally, the numerical-integration 
formulas resemble those of the Lax-Wendroff method. A com prehensive analysis 
of our numerical technique was given by Kshevetskii ( 2006f l. 

Our simulations are carried out under the assumption that the velocity com¬ 
ponent normal to the boundary of the computation domain and the normal 
derivative of the magnetic field vanish at this boundary. The computation do¬ 
main is chosen to be large enough to eliminate boundary effects affecting the 
processes in its central part, where the phenomena under study develop. 

We solve the problem for a chromospheric layer where the temperature of both 
the ion (proton) and electron components of the plasma is 25,000 K, which gives 
a total characteristic temperature of T* = 50, 000 K f or the plasma medium. 
We assume here W = 10^® m“^ (Demoulin and Klein, l200r)h . In such a layer, 
0* = 1.2 X 10“®, K* = 3.8 X 10“"^, and = 5 x 10“® (see Appendix A). 

Further, we consider the case of /3o* = 1.6 turning sometimes to those of 
/3o» = 1-5 or ,00* = 2.3. In accordance with Section |31 the dimensional units for 
the case of ,5o* = 1-6 are: Bq = 0.33 G, 


vo» = 23 kms 


to* = 1.1 minutes. 


(17) 


The applicability of the our initial-value problem to description of plasma 
dynamics at this intensity of Bq is discussed in Appendix D. 
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The calculations for M = 0 were carried out in a domain with Xmax = 0.207, 
Zmax = 2.8, which corresponds to 300 km (320 grid points) in height and 4200 km 
(988 grid points) in the horizontal direction. For M ^ 0, the domain measures 
0.207 X 0.639, or 300 x 922 km^, with 80 x 247 grid points. 


6. Results 


6.1. The Phenomenon of Velocity Burst 


Our simulations reveal an ensemble of primary waves produced by contraction 
and expansion of the whole current zone associated with the unbalanced com¬ 
ponent of B (see Section and Figures [T] d,c). Against the background of these 
waves, at the centre of the magnetic configuration, the pinch-sausage instability 
begins developing and becomes especially pronounced after the waves leave this 
area (see Section [7]). 


As shown previously (Alekseeva and Kshevetskii. l2011r ). the nonlinear stage 
of the pinch-sausage instability forms cross-shaped transient structures, which 
are more complex and smaller in size in the calmer PIB case, as compared with 
the more dynamical, completely unbalanced NIB case. The development of the 
pinch-sausage instability in the PIB case passes through two early stages: only 
after this pinch effect has manifested itself in the unbalanced component of 
the magnetic field and has produced its cross-shaped pattern, does the sausage 
instability set in in the balanced component of the field, and the concentration of 
the whole magnetic field approaches the symmetry axis x = Xc of the initial mag¬ 
netic configuration. This gives rise to intricate transient flow structures, which 
facilitate the development of turbulence and, accordingly, favour the heating 
of the plasma. Thus, T can be observed to double in the PIB case, while no 
appreciable temperature changes are associated with the pinch instability in the 
NIB case. 

Now, simulations based on the improved algorithm ensuring higher accuracy 
reveal another peculiarity of the coevolution of the magnetic field and plasma 
state in the PIB case (Figure^]). Against the background of the above-mentioned 
fairly gradual variations in the physical quantities, impulsive phenomena of 
increasing plasma velocity, i.e. velocity burst, develop, no similar phenomena 
being observed in the case of NIB. 

The numerical result obtained here agrees with the properties of the exact 
solution deducible from the mathematical form of our initial-value problem (see 
Section!?]). Thus, in the PIB case, we observe the whole scenario of the sausage 
pinch instability, up to the manifestation of the final velocity burst. (We note 
that no exceptions occur in the process of computation, and the code continues 
running with progressively decreasing steps). 

Generally, /3o* is higher, the slower the process (Figure |2|); we also observe 
this regularity in other numerical simulations, not described here, although there 
is no smooth dependence of the parameter values /Iq* (see Figures [7^, b for an 
example). 
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Figure 2. The maximum absolute velocity over the computation domain as a function of time 
in the cases of PIB (solid curve) and NIB (dotted curve) at ^o* = 1-5 (a), /3o* = 1-6 (b) and 
/5o* = 2.3 (c) in the absence of gas viscosity [M]. 
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Physical reasons for this effect can naturally be expected, since we solve 
an initial-value MHD problem with an unstable initial magnetic configuration, 
where the sausage instability develops under different physical conditions due to 
the presence of interfering primary waves, even at very close values of /3o* (see 
Section 17.fl below). 

Also, no smooth transition between the final stage of the numerical solutions 
should be expected since the required accuracy can more hardly be achieved. 
For example, in both cases /3o» = 1-5 and 1.6, our numerical solutions represent 
the evolution of the process with the gradual velocity variations changing into a 
dramatic increase (the white areas in Figuresl^ti-jb). This increase should result in 
a rapid growth of both physical and numerical instabilities with the spatial scale 
of disturbances approaching zero. At this stage (the shaded areas in Figures[2^,b 
right of about t = 0.27), the simulation requires progressively shorter time 
steps and catastrophically loses its accuracy. Thus, the white and shaded areas 
correspond to a normal and a poor accuracy, respectively (this is why the small 
change in /3o* can delay the final infinite velocity growth). However, we do show 
the curves in the shadowed areas in Figures [2^, b, since they demonstrate that 
our numerical solution actually reproduces the analytically predicted velocity 
divergence. In no way does it give a quantitative description of the final stage of 
the processes. 


6.2. Velocity Burst and the Known Pinch-Effect Features 


We consider the process in more detail precisely for the case /3o* = 1.6, where 
individual peaks can be singled out. Each of them has features of an underde¬ 
veloped final velocity burst. The solution at /3o* = 1.6 shows that a first, faint 
velocity burst with a peak velocity of vi = 1.2 at t = 0.224 precedes the second 
one with V 2 = 3.9 at t = 0.266 (Figure [2|). The first burst appears when the 
cross-shaped pattern has already formed (Figure [3]). 

Later, the distributed magnetic field approaches the symmetry axis x = Xc 
to occupy an area initially free of magnetic field (Figure 01 left). This induces 
a sharp increase in the current-aligned velocity (which is seen as a brightening 
along the axis x = Xc in Figure 01 right). The high-speed jets corresponding to 
the peak of the velocity burst are located in the central part of the magnetic 
configuration (Figure[51). This velocity burst occurs in a time interval of duration 
0.0028. 


Since the plasm a velocity v in the system of MHD equations (Brushlinskii 
and Morozov, 1980h coincides with the bulk velocity of the ion component Uj, at 
the velocity burst, an ion acquires an additional kinetic energy of jl. This 
corresponds to the ion kinetic temperature T^i = 511,000 K gained rapidly, in 
0.1 seconds. Thus, we see that the velocity burst produces suprathermal protons. 


In Section O we assumed that the gravitational acceleration of plasma 
should not be taken in to account in our consideration. Now we can substantiate 
this assumption knowing that large velocities are achieved (Figure [21). Even in 
the extreme case of free fall, in a time interval of At ^ 5.3 seconds, a plasma 
parcel would gain a velocity increment of Au ^ glSi ^ 1.45 kms“^, which 
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Figure 3. Distribution of B (left) and v (right) in the plane (x,z) at t = 0.216, M = 0. 



Figure 4. Distribution of S at t = 0.256 (left) and v at the onset of the velocity burst, 
t = 0.264, (right) in the plane {x,z). 

is vanishingly small compared to the typical plasma velocities, not to mention 
the burst velocities. (Including this estimate was kindly suggested to us by the 
referee.) 

6.3. The Increase of the Magnetic-Pressure Gradient as the 
Producer of the Velocity Burst 

The idea of the physical mechanism of the velocity burst can be formed if we 
consider the behaviour of other quantities during this burst. Only the region 
shown in Figure [SJ right, will be considered in what follows. The subscripts max 
and min will refer to maximum and minimum values of a quantity at hand over 
this very region. 

Note that B, p, V and T vary in the burst area not so strongly as v does (Fig¬ 
ure (5]). However, according to Equation ([5]), the velocity variations are directly 
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Figure 5. Distribution of v (left) in the plane (x^z) at the time of the second velocity peak, 
t = 0.2656. Right: vector map of v (magnified) for the same time. A division of the vertical 
(a:) axis is 0.93 km; of the horizontal (z) axis, 4.25 km. The longest arrow corresponds to a 
velocity of V 2 = 3.9. 



t 


Figure 6. Temporal variation of MHD quantities near the time of velocity burst: Umax (heavy 
solid curve), (-B^/2)max (light solid curve), 2.5Tmax (heavy dotted curve), 2.5Tmin (light 
dotted curve), 2.5pmax (heavy dashed curve), 2.5pmin (light dashed curve), 2.57^max (heavy 
dot—dashed curve), 2.5'Pmin (light dot-dashed curve). 


controlled by the gradient of the fnll pressure rather than the above-mentioned 
quantities. Note that the elevated-gradient areas of the magnetic and the full 
pressure are virtually coincident (Figure O, the former giving a dominant con¬ 
tribution to the latter. Consider how the magnetic-pressure gradient is related to 
the velocity (FigureO- Indeed, the velocity vector is aligned with the gradient of 
the magnetic field pressure wherever both are graphically representable, although 
this is not seen everywhere. However, the gradients cannot always be judged by 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 17 




















L.M. Alekseeva, S.P. Kshevetskii 






5 


5 

♦ 

^ 0 

# - 

^ 0 

4 #• 

❖ 

-5 


-5 



-20 -10 0 10 20 -20 -10 0 10 20 

z z 


Figure 7. Maps of the total pressure (left) and magnetic pressure (right) at an intermediate 
development stage of the burst (f = 0.2648). Areas where the pressure exceeds its mean value 
over the region are hatched. 


the contours plotted; this is all the more so because even large gradients can be 
handled by our code using generalised functions. 

The velocity-burst mechanism could likely be described as follows. In such an 
impulsive process, the velocity increase achieved over a time interval [dt] is 


8v 


St- 




(18) 


Generally, |Vi3^/2| can reach large values, and our technique makes it possible 
to simulate cases with large gradients. 

A link between velocity increases and the presence of large magnetic-field 
gradients can be revealed by an analysis of temperature variations in the burst 
area. The temperature is increasing during the burst (Figure [5]); this occurs 
where the density values are reduced and the temperature increase is therefore 
not an adiabatic effect of gas compression (Figure [9|). Actually, the plasma is 
heated by electric currents, and the temperature is the higher, the larger the 
current density. In our 2D geometry, the current-density magnitude is propor¬ 
tional to |Vi3|. A comparison of the temperature-velocity maps (Figure [10]) for 
the initial stage of the process and for a well-developed burst shows that the 
longest arrows issue from the areas of increased temperature, i.e. of increased 
gradients of B. 

The velocity burst revealed here is accompanied by large velocity gradients; 
therefore, this burst should be affected by viscous forces if M 0 is taken 
into account in Equation ([S]). To verify our suggestion that the velocity-burst 
mechanism stems from Equation (1181) . we have carried out calculations for the 
case of M 7 ^ 0 assuming /i* value exaggerated by a factor of 20 compared with the 
actual viscosity calculated according to Equation dMj) for the considered solar- 
plasma layer. The presence of such a viscosity was even able to suppress the 
dramatic growth of Umax at the final development stage of the pinch instability. 
The evaluation of the specific peak value of Umax requires taking into account 
the dependence of viscosity on the local temperature and using an accurate 
expression for the viscous-stress tensor. 
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Figure 8. Normalised velocity ujv and contours of the magnetic pressure for an early 
(t = 0.2640, left) and the highest (t = 0.2656, right) development stages of the burst. The 
dotted curve corresponds to a level of (l/3)[(B^/2)max — (-S^/2)min]- The unit-length arrow is 
shown near the lower left corner of each panel; it corresponds to the longest arrow in Figure [5] 
right. 



Figure 9. Normalised velocity vn, temperature T(solid curves) and density p (dotted curves) 
for the highest stage of burst development [t = 0.2656). Contour levels from 1.456 to the peak 
value with an increment of 0.076 are used for T; from 0.736 with an increment of 0.086, for p. 
The unit-length arrow is shown near the lower left corner; it corresponds to the longest arrow 
in Figure [S] right. It can be seen that the temperature reaches its maximum where the density 
is reduced. 


7. Discussion 

7.1. The Onset of the Sausage Instability Against the Background of 
Interfering Waves 

Thus, our numerical simulations reveal a velocity burst, which is observed only 
at the initial conditions of the PIB type, the case physically corresponding 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 19 



































































































L.M. Alekseeva, S.P. Kshevetskii 




Figure 10. Normalised velocity vn and temperature for an early (t = 0.2640, left) and the 
highest (t = 0.2656, right) development stages of the burst. The lowest contour level 1.28 and 
the contour increment 0.048 are the same in both panels. The unit-length arrow is shown near 
the lower left corner of each panel; it corresponds to the longest arrow in Figure [5] right. 



Figure 11. The maximum absolute velocity over the computation domain as a function of 
time in the case of PIB at different viscosity /i* = 0 (solid curve) and /r* = 10“^ (dashed 
curve). 


to the penetration of a new, unbalanced magnetic field to the area of pre- 
existing field, which has already reached pressure balance ISection I3.2I1 . This 
impulsive phenomenon is given by the spontaneously developing sausage in¬ 
stability in the course of the general plasma and magnetic-field evolution (see 
Section 113). At some time, the instability begins developing in the central part 
of the configuration amidst the primary waves produced by the unbalanced B 
(Figure [12]). Ultimately, it leads to the divergence of the solution, which was 
predicted analytically. 

The rate of the instability development should depend on the spatial dis¬ 
tribution of the quantities at the onset of instability (as the growth rate of the 
amplitude depends on the perturbation wavevector in equilibrium configurations 
according to the linear theory of the sausage instability). Therefore, the partic¬ 
ular interference pattern of primary waves at the time of onset should affect 
the particular scenario of instability development and, eventually, the time of 
divergence. 
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Figure 12. The maximum value of the velocity component along the contact zone over the 
computation domain in the case of PIB /3o* = 1-5 at £ = 0.08 (left) and £ = 0.136 (right). The 
sausage instability starts developing amidst primary waves. 


This obvious feature of the solutions should be kept in mind comparing the 
solutions for close values of the external parameters of the problem, e.g., for 
/?o* = 1-5 and 1.6 (Figures [2^,b). Although processes of the type ii) giving 
rise to primary waves (see Section IQ) differ very little, stronger magnetic 
fields ensure greater amplitudes and steeper profiles of waves; therefore, the 
wavevectors of the interference patterns is highly sensitive to the parameter /3o* 
variations. This results in different temporal variations of the maximum velocity 
seen in the white regions of Figures [2^ and[2jD (note that, at different times, the 
maximum velocity is generally reached at different points, directly depending 
on the temporal variation of the interference pattern). At later times, when the 
disturbance related to the sausage instability becomes strong, the differences in 
the behaviour of the curves become greater (the shaded regions of Figures [5^ 
andj^Jo). Thus, the manifestations of the velocity burst depend on the fine details 
of the early development stage of the instability. 

7.2. Velocity Burst and the Problem of Dynamic Chromospheric 
Energisation 

In our initial-value problem, the plasma is initially motionless and has a constant 
temperature of 50,000 K throughout the plasma volume. Further, the magnetic- 
field energy transforms into the thermal and kinetic energy of the plasma because 
of the development of a pinch with a sausage instability. 

For example, in the numerical simulations at /3o* = 1.6 in the absence of 
gas viscosity, the plasma temperature increases by a factor of two, whereas the 
kinetic temperature of the proton component reaches a peak value of 0.511 MK 
at undimensional t = 0.2656, or an average particle energy of 44 eV, and sub¬ 
stantially exceeds it at later times during the final velocity burst (Figure [2)o; 
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we note that the exact solution goes to infinity-see Section n. We see that 
impulsive processes can manifest itself against the background of a gradual 
development of the pattern. Specifically, intense plasma streams emerge in a 
temporal interval of about 0.1 second or less near the centre of the magnetic 
configuration. Such a velocity burst can physically be interpreted as an event of 
suprathermal-ion generation. This is indicative of the dynamic energisation of a 
parcel of chromospheric plasma to coronal temperatures and, if the final stage 
is achieved, of the increase in the kinetic temperature to high values, even those 
observed in flares. 

Note that magnetic reconnection is not present in our solution; it would 
generate plasma flows varying along the magnetic-field line, which is forbidden 
by the chosen 2D geometry without variations of the physical quantities in the 
direction of B. Therefore, under the conditions of the upper chromosphere, the 
pinch-sausage effect can be considered a possible alternative to the magnetic- 
reconnection process as the producer of flares. However, it can develop only in 
the presence of the “old” magnetic field, which has already reached pressure 
balance (we recall that, in Figure O the curves for PIB case, in contrast to NIB 
case, demonstrate the velocity burst; see also Appendix C). 

To describe the velocity burst more accurately, the equation system used 
here should be generalised by adding more physical processes. Thus, our simu¬ 
lations suggest that a finite gas viscosity can reduce the kinetic temperature of 
the fast streams. There are a number of possibilities to obtain more realistic 
solutions: introducing an adequate gas viscosity in the system of collisional 
MHD equations; including terms with the Hall-plasma-dynamics parameter 
i.e. considering a conductivity tensor with cth and crp present; and proceeding 
to three-dimensional simulations (which, in particular, could make it possible to 
consider also magnetic reconnection). This work is currently underway. 


7.3. Velocity Burst and the Ion Beams 


The suprathermal ions (protons) revealed in our simulations are related to a phe¬ 
nomenon known from laboratory nuclear-fusion experiments probably. Specifi¬ 
cally, neutrons were produced in Z-pinch plasma by relatively small quantities 
of “beam” deuterium ions accelerated in the direction of the current to energies 
of 50 - 200 keV due to the development of the “sausage” m = 0 instability of the 
pinch. Colliding with the deuterium ions whose temperature was much lower, the 
beam ions produc ed fus ion neutrons, which were therefore not thermonuclear 
(Velikovich et al., 2003)- The origin of this phenomenon due to the nonlinear 
evolution of the m = 0 instability of a stagnated compressional Z-pinch is still not 
fully understood. Some of the proposed mechanisms are fluid-like and others are 
kinetic in nature. From the MHD standpoint, ion jets in both directions should 
occur equally from the “necked” region. Presumably, it requires introduction of 
the Hall term and/or, the finite ion Larmor radius terms in the stress tensor to 
break th is sym metry and allow a beam in the direction of the current to form 
(Haines, l201lh . 

Some results of our numerical experiment with MHD hydrogen proton plasma 
show evidence of their correspondence to this situation. In the PIB case, where 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 22 








MHD Simulation of Chromospheric Plasma 


stagnation is initially present, a nonlinear state of the pinch process gives rise to 
opposite current-aligned ion jets of short duration in the centre of the magnetic 
configuration (Figures [5] and [5]). The energies of ions accelerated during the final 
velocity burst would be comparable with the energy value indicated above for the 
deuterium ions obtained in laboratories. It is worth adding that the simulated 
pair of suprathermal-ion jets could actually produce an ion beam accelerated in 
the direction of the current. Indeed, suprathermal ions experience less frequent 
collisions than the thermal particles of the plasma and are more directly con¬ 
trolled by the electromagnetic held. Near the symmetry axis along which they 
move, the magnetic held is weak. The electric held should accelerate the ions 
of one jet and decelerate the ions of the oppositely directed jet. Therefore, only 
one ion beam will be produced. 


8. Conclusion 

Contacted oppositely directed horizontal magnetic helds are not rare in the up¬ 
per chromosphere, where the plasma is fully ionised. Therefore, the phenomena 
described by our numerical solutions should not be unusual. 

We have traced the coevolution of the plasma and magnetic held from the 
origin of the steadying waves to the hnal stage when a velocity burst occurs. 
The temperature of the chromospheric plasma increases by two to three times in 
the course of the evolution. The velocity burst is due to the sausage instability 
that develops against the background of the steadying waves, turbulence, frontal 
zones, etc. As a result, fast streams suddenly appear at the centre of the magnetic 
conhguration. The kinetic temperature corresponding to such fast streams can 
approach the values recorded in the corona and, likely, in hares. So, the pinch- 
sausage effect can be a possible energiser of the upper chromosphere and an 
alternative to the magnetic-reconnection process in producing hares. 

We have numerically solved the standard self-consistent MHD system of equa¬ 
tions. This made it possible to relate the plasma-magnetic-held coevolution to 
the initial values of the basic MHD quantities without invoking ad-hoc non¬ 
observable parameters. The standard form of the MHD system allows us to 
compare our results for the solar plasma with those for the laboratory plasma, 
and we note that the velocity burst revealed in our simulations resembles the 
sudden emergence of high-speed proton beams in laboratory Z-pinch devices. 


Appendix 


A. Material Properties of the Medium 


Magnetic diffusivity. 
conductivity [cr] of a 


Following Spitzer (19621) and Priest ( 1982 1. we dehne the 
collisional, fully ionised hydrogen plasma as 


e^NeTe 

me 


(19) 
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where Ng and Tg are the local electron concentration and electron-ion collision 
time. In a fully ionised gas there is so me unc ertainty as to the appropriate value 
of Te to use in this equation (Spit zer, Il962h . We choose in accordance with 
formu la (2.5e) of Braginskii ( 1965ll . which coincides with Te in § 5.7 of B alescu 
1 1988f l. [Thus, ou r a is equal to that used by Brushlinskii and Morozov ( IQSOl l 
and Brushlinsky (198^. It is also used as a basis for the generalisation of the 
form o f con ductivity to the case of plasma with neutral particles (Leake and 
Arber, 2006ll .] Explicitly writing Te yields 0 in the form Equation (IIOL where 


TOgWi c 


1 Ae^ 


27r ff {kT^f 0.75 


( 20 ) 


(A is the Coulomb logarithm, which we set here equal to 20). 

In a quasi-neutral plasma, a oc NeTe is a function of only the local tempera¬ 
ture. 

Thermal conductivity and gas-dynamic viscosity coefficient. 

The thermal conductivity of a collisional plasma differs substantially between 
the region of magnetised plasma and the region o f zero magnetic field. In the 
(dimensional) heat-transfer equation (Braginskii, Il965l l. we use the quantity 
corresponding to the case of weakly magnetised plasma (where, in the com¬ 
monly accepted notation, k± ^ K(j_ ~ k|). Once the nondimensionalisation of 
the equations is accomplished, this coefhcient, appearing in the expression for 
Equation determines 


_ ksT^'^y/rni 

~ 32i7fc3/2A^* ’ 
3.16-0.75 
Av^ ^/^ ■ ’ 


( 21 ) 

( 22 ) 


(note that the so-called Spitzer thermal conductivity [fcg] has dimensions of 
ergs“^ cm“3 deg”^'^^). This substantially overstates the thermal conductivity 


in the region of the str ong m agnetic field; however, as our results indicate (Alek¬ 
seeva and Kshevetskii. I 2 OI 111 , major thermal changes occur in the region of weak 
magnetic field x = Xc- 

The expression for the kinematic-viscos ity co efficient of plasma with a mag¬ 
netic field is very complicated (Brag inskii, Il965l l. Instea d, even in the case of a 
magnetic held present. Priest (|l982[ l and Spitzer (119621 ) consider the coefhcient 
of plasma without magnetic held 


Mo = 2.21 X 10"^^A-^(ri/K)^/2g • cm 




(23) 


where T{ is the ion temperature in K (since the gas-dynamic viscosity of plasma 
is determined by ions, it is worth recalling that the ion temperature is half the 
plasma temperature under consideration). That practice is appropriate for use 
in our investigation, since, in the magnetic helds studied with the pinch effect, 
the most interesting phenomena arise in the neighbourhood of the neutral line 
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where magnetic field is not great. So we obtain M in Equation (USD with the 
dimensionless parameter [/i*] determined by 


0.215 T,/(l K) 

A TV* X (1 cm^) 


(24) 


The present study addresses the qualitative role of viscosity; therefore, the 
simplest form of the viscous force is used in Equation ([5]). 


B. Initial Magnetic-Field Distribution 

The structure of the magnetic held at the initial time is determined by the 
expression 

B{x,z)\t=o = D{x,z) + d{x,z), (25) 

D = [l + Cexp(-C)]tanh[g(a;-a;c)], C, = {z - z^f jzi, (26) 
^ ^ ^ s(a;-Xc)^exp(-1000C) 

[xl + {x-x,Yf[l + s,{z-z,Yy 

where ([ — 130, Xg — 0.031, Xc — ^max/2; — ^niax/2. 

Plasma with zero gas-dynamic viscosity was simulated in a computational 
domain with Xmax = 0.207, z^ax = 2.8 at si = 1.76, Zk = 0.17, s = 0.01. 
With this choice of constants, the term |I1| in Equation (1251) is approximately 
half of the value of |(i|. For a viscous plasma, the calculations were carried out 
in a computational domain with Xmax = 0.207 and Zmax = 0.639 at si = 0, 
Zfc = 0.025, s = 0.00025. Therefore, we get |I1| ~ |d|. 


C. Initial Pressure and Density in the PIB Case 

Recall that p(x, z)|t=o = 1 outside the region occupied by the magnetic held. As¬ 
sume that, in the region of the magnetic held Equation (1251) . only the component 
£> is in equilibrioum with the gas pressure: 

r(x,z)lt=o + D(x,zf/2 = /3o,/2. (28) 

Taking into account Equation © leads to the following initial conditions for 
pressure and density: 


V(x,z)jt=o = /3o*t2- I)(x,zf/2, 

(29) 

1 1 D(x,zf 

'’I- = ‘ /Jo. ^ 

(30) 
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D. The Applicability of the Considered Initial-Value Problem 
to the Case of /3o* = l-S — 2.3 

Weak magnetic field evolves slowly and requires computations over long periods 
to reach the divergence of the solution (described in Section |4|) . This leads to 
accumulation of errors. Therefore, we choose the characteristic B value to be as 
large as Bq = 0.33 G. 

Although this Bq corresponds to characteristic values uJeoTe* 1, we nev¬ 
ertheless use it to solve numerically the initial-value problem - Equations © 
- ®, (HU), (1^ and either (fl^ . (ITI)ll or (1^ . (fliTll - with isotropic 
intending to supplement the system, at a later time, with terms including the 
parameter C or, in other words, to automatically take into account the anisotropy 
of the local coefficients in the areas where this anisotropy is substantial (see 
footnote in Section El). 

However, in the magnetic configuration considered (Figure [1^), the most im¬ 
portant phenomena develop near the line of contact of opposite polarities x = Xc, 
where magnetic field vanishes and the local coefficients are therefore isotropic 
(see SectionEl). Thus, the results obtained for the given Bq appear to adequately 
reflect the reality. (We can also note that our tentative calculations with ^ 0 

for the same magnetic configuration and initial conditions show virtualally no 
difference with the case of ^ = 0 presented here; only a more rapid evolution is 
observed.) 

Acknowledgments We cordially thank A.V. Getling for carefully considering the ideas of 
this study and useful recommendations. We are also grateful to the referee for discussions. This 
work was supported by the Russian Foundation for Basic Research (project no. 12-02-00792-a). 


Disclosure of Potential Conflicts of Interest 

The authors declare that they have no conflicts of interest. 


References 

Alekseeva, L.M.: 1980, Physical nature of the anode explosion in plasma channels. Sov. Tech. 
Phys. Lett. 6(11). 

Alekseeva, L.M.: 1997a, Generation of magnetic field by acceleration of Hall plasma and some 
some consequences of this process. In: Bordage, M.C., Gleizes, A. (eds.) XXIII Internat. 
Conf. on Phenomena in Ionized Gases 2, Universite Paul Sabatier, Toulouse, 156. 

Alekseeva, L.M.: 1997b, Hall effect as a producer of effective gravity aand of quasi- acoustic- 
gravity waves in plasma flows. In: Bordage, M.G., Gleizes, A. (eds.) XXIII Internat. Conf. 
on Phenomena in Ionized Gases 1, Universite Paul Sabatier, Toulouse, 206. 

Alekseeva, L.M.: 1999, Effective gravity due to the Hall effect and quasi-acoustic-gravity waves 
in accelerated plasma flows. J. Plasma Phys. 61, 671. DOI ADS 

Alekseeva, L.M.: 2006, The Hall effect as the producer of current sheets in astrophysical 
plasmas. Astron. Lett. 32, 274. DOI ADS 

Alekseeva, L.M., Kshevetskii, S.P.: 2006, Numerical simulations of the MHD structuring of 
rarefied solar plasma. In: Multiwavelength Investigations of the Sun and Modern Problems 
of Solar Activity, All-Russia Gonference, Abstracts, Special Astrophysical Observatory, 
Russian Academy Sci., 12. 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 26 



MHD Simulation of Chromospheric Plasma 


Alekseeva, L.M., Kshevetskii, S.P.: 2011, Solar plasma structuring due to two-dimensional 
pinch effect over the photosphere. Astrophys. Bull. 66(1), 72. 

Alekseeva, L.M., Solov’ev, L.S.: 1964, Current votices and critical surfaces in a magnetihydro- 
dynamic flow. Prikladnaya Matematika i Mekhanika 28(6), 987. 

Aschwanden, M.J.: 2001, An evaluation of coronal heating models for active regions based on 
Yohkoh, SOHO, and TRACE observations. Astrophys. J. 560, 1035. DOI ADS 
Aschwanden, M.J.: 2004, Physics of the Solar Corona. An Introduction, Springer, Berlin. 
Aschwanden, M.J.: 2008, An observational test that disproves coronal nanoflare heating 
models. Astrophys. J. Lett. 672, L135. DOI ADS 
Aschwanden, M.J., Nightingale, R.W.: 2005, Elementary loop structures in the solar corona 
analyzed from TRACE triple-filter images. Astrophys. J. 633, 499. DOI ADS 
Balescu, R.: 1988, Transport Processes in Plasmas 1, Elsevier, Amsterdam. 

Braginskii, S.I.: 1965, Transport processes in a plasma. In: Leontovich, M.A. (ed.) Rev. Plasma 
Phys. 1, Consultants Bureau, New York, 205. 

Brushlinskii, K.V., Morozov, A.I.: 1980, Calculation of two-dimensional plasma flows in chan¬ 
nels. In: Leontovich, M.A. (ed.) Rev. Plasma Phys. 8, Consultants Bureau, New York, 105. 
ADS 

Brushlinskii, K.V., Gerlakh, N.I., Morozov, A.I.: 1968, Two-dimensional steady flow of highly 
conducting plasma in a coaxial system. Fluid Dynamics 1, 134. DOI ADS 
Brushlinsky, K.V.: 1989, Computational models in plasma dynamics. In: Dwoyer, D.L., Hus- 
saini, M.Y., Voigt, R.G. (eds.) Numerical Methods in Fluid Dynamics, Lecture Notes in 
Physics 323, Springer, Berlin, 21. ADS 

Carlsson, M.: 2007, Modeling the solar chromosphere. In: Heinzel, P., Dorotovic, I., Rutten, 
R.J. (eds.) The Physics of Chromospheric Plasmas CS-368, Astronomical Society of the 
Pacific, San Francisco, 49. ADS 

Demoulin, P., Klein, K.-L.: 2000, Structuring of the solar plasma by the magnetic field. 
In: Rozelot, J.P., Klein, L., Vial, J.-C. (eds.) Transport and Energy Conversion in the 
Heliosphere, Lecture Notes in Physics 553, Springer, Berlin, 99. ADS 
Fletcher, L., Turkman!, R., Hudson, H.S., Hawley, S.L., Kowalski, A., Berlicki, A., Heinzel, P.: 

2010, Solar flares and the chromosphere. ArXiv e-prints: 1011.4650v2 [astro-ph.SR]. ADS 
Gabriel, A.H.: 1994, Some problems in understanding the solar corona. In: Rusin, V., Heinzel, 
P., Vial, J.-C. (eds.) lAU Colloq. Iff: Solar Coronal Structures, Slovak Academy of 
Sciences, Tatranska Lomnica, 1. ADS 

Gary, G.A.: 2001, Plasma beta above a solar active region: Rethinking the paradigm. Solar 
Phys. 203, 71. DOI ADS 

Goodman, M.L.: 2005, Self-consistent magnetohydrodynamic modeling of current sheet struc¬ 
ture and heating using realistic descriptions of transport processes. Astrophys. J. 632, 1168. 
DOI ADS 

Haines, M.G.: 2011, A review of the dense Z-pinch. Plasma Physics and Controlled Fusion 
53(9), 093001. DOI ADS 

Kshevetskii, S.P.: 2006, Study of vortex breakdown in a stratified fluid. Comput. Math. Math. 
Phys. 46, 1988. DOI ADS 

Kshevetskii, S.P., Gavrilov, N.M.: 2005, Vertical propagation, breaking and effects of nonlinear 
gravity waves in the atmosphere. J. Atmos. Solar-Terr. Phys. 67, 1014. DOI ADS 
Leake, J.E., Arber, T.D.: 2006, The emergence of magnetic flux through a partially ionized 
solar atmosphere. Astron. Astrophys. 450, 805. DOI ADS 
LeVeque, R.J.: 1992, Numerical Methods for Conservation Laws, 2nd edn. Birkhauser, Basel. 
Morozov, A.L, Solov’ev, L.S.: 1980, Steady-state plasma flow in a magnetic field. In: 

Leontovich, M.A. (ed.) Rev. Plasma Phys. 8, Gonsultants Bureau, New York, 1. ADS 
Priest, E.R.: 1982, Solar Magnetohydrodynamics, Geophysics and Astrophysics Monographs 
21. ADS 

Schrijver, G.J.: 2001, The coronae of the Sun and solar-type stars (GD-ROM Directory: 
contribs/schrijv). In: Garcia Lopez, R.J., Rebolo, R., Zapaterio Osorio, M.R. (eds.) 11th 
Cambridge Workshop on Cool Stars, Stellar Systems and the Sun CS-223, Astronomical 
Society of the Pacific, San Francisco, 131. ADS 
Spitzer, L.: 1962, Physics of Fully Ionized Gases, Interscience, New York. 

Thomas, J.P., Roe, P.L.: 1993, Development of non-dissipative numerical schemes for computa¬ 
tional aeroacoustics. In: AIAA 11th Computational Fluid Dynamics Conference, Technical 
Papers 2, Orlando, FL, 906. ADS 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 27 



L.M. Alekseeva, S.P. Kshevetskii 


Velikovich, A.L., Clark, R.W., Davis, J., Chong, Y.K., Deeney, C., Coverdale, C.A., Ruiz, 
C.L., Cooper, G.W., Nelson, A.J., Franklin, J., Rudakov, L.L: 2007, Z-pinch plasma neutron 
sources. Phys. Plasmas 14(2), 022701. DOI ADS 


SOLA: Alekseeva.tex; 17 August 2015; 0:24; p. 28 



